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Abstract — In the present study, application of Amberlite FPA-9& was investigated for the removal of Acide 
Orange 7 from aqueous solution using I lie continuous method and was optimized using Box-Behnken design 
(BBD) and full factorial design (FFD). Fixed bed adsorption has'become a frequently used in wastewater 
treatment processes, ffi^this work, the intention of the study was tp explore the efficacy and feasibility for azo 
dye, Acid Orange 7 (A07) adsorption onto fixed bett~columimf Amberliti* FPA-98. The effect of operating 
parameters such as flow rate, fiuititu dye concentration, and bef height jvas modeled _ by response surface 
methodology (RSM). This’tiucJy. compares Box-Behnken design (BBD) and full factoriaT design (FFD) utility for 
modeling and optimization by''response surface methodology. The precision of the equation obtained by RSM 
was confirmed b\j the analysis bf variance (ANOVA)and calculation of correlation coefficient relating the 
predicted and the ffperimentctl values of adsorption efficiency. The res^LFtrc^vealed a good agreement between 
the predicted values, as obtained by full factorial design (FFD) a id the experimental values for A07 (R2 = 
0.987) . The optimum conditions proposed by Box-Behnken design (BBDf'k) reach the maximum dye removal 
through adsorption process. Under the optimum conditions, the removal efficiency of A07 were (R2 = 0.959). 
The application of respohseiyurface methd4ology in order to optimize using Box-Behnken design (BBD) and full 
factorial design (FFD). The tTsedVch- On moaFl\ng adsorption by RSM has been highly developed and The 
Amberlite was shown to be suitable adsorbent for.adsorption ofA07 using fixed-bed adsorption column. 
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I. INTRODUCTION 

T extile industry is very greedy -in i .wat’ef--and r tITus~[b] gcttdralCs an impQEtaTit-jquantity of effluents highly 
charged with pollutants which constitute "d serious threat for the environment. Consequently, these effluents 
require a preliminary treatment in order to decrease their polluting load before being rejected into the natural 
environment. It is considered that the textile i ndustry js r.eSponsiblc for 15% to 20% of the global water pollution 
[2]. Among the discharged pollutants, organic dyes'of are n<3t Only rtTsponsiblcfor an esthetic pollution of water, 
but also count among the most toxic compounds, even at low concentration [3]. Today, about 10,000 different 
dyes are produced worldwide, for a global production of 7.105 tons per year. 

Dye effluents discharged from the dyestuff manufacturing, dyeing, printing and textile industries 
represent a serious problem all over the world. They contain different types of synthetic dyes which are known 
to be a major source of environmental pollution in terms of both the volume of dye discharged and the effluent 
composition [4]. Most of these dyes are toxic, mutagenic and carcinogenic [5]. From an environmental point of 
view, the removal of synthetic dyes is of great concern. 

Natural colorants have been used since prehistoric times. In 1856, Perkin's discovery of mauve marked the start 
of the modern synthetic dye industry. In the last 145 years, several million different colored compounds have 
been synthesized, with ca. 15,000 colorants over time, produced on a commercial scale. The annual worldwide 
production of dyes is approximated at 800,000 tonnes and about 50% of these are azo dyes [6]. Anionic azo dyes 
contain many compounds from the most varied classes of dyes, which exhibit characteristic differences in 
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structure (e.g. azoic, antraquinone, triphenylmethane, and nitro dyes) possessing water-solubilizing, ionic 
substituents as a common feature. Acid Orange 7 has been reported to induce bladder tumors. It can also easily 
undergo enzymatic breakdown along with reduction and cleavage to give aromatic amines, which, upon 
exposure, can cause methemoglobinemia. When Acid Orange 7 enters the human body through ingestion, it is 
considered genotoxic; however, if some impurities, such as aromatic amines, are present, it shows mutagenic 
activity. Due to large amounts of Acid Orange 7 consumption, it is essential to have a proper method to remove 
this dye from wastewater in order to avoid potential threat for the environment. As synthetic dyes in wastewater 
cannot be efficiently decolorized by traditional methods, also Acid Orange 7 does not decompose biologically, 
and resists to light irradiation and chemical oxidation. 

In general, chemical and physical discoloration methods have been used for dye removal from 
wastewaters [7] such as coagulation and flocculation processes which are largely used for wastewater treatment 
in the textile industries. However application of these methods is somewhat restricted due to some limitations 
such as operational costs, formation of hazardous by-products, intensive energy requirement and limited 
adaptability to a wide range of effluents and these processes are not always effective for dye removal. Moreover, 
they can sometimes generate secondary pollution due to the excessive use of chemical reagents. 

However, the adsorption of this dye on efficient solid supports is considered as a simple and economical method 
for its removal from water and wastewater providing sludge-free cleaning operations and many studies have 
been conducted to find suitable adsorbents to reduce Acid Orange 7 concentration [8]. 

Because a number of exchange resins have been used quite efficiently for the removal of specific organic 
compounds [9], this study investigates the adsorption characteristics of Acid Orange 7 dye on the strongly basic 
polyacrylic anion exchanger Amberlite FPA-98 of macroporous structure. The present work deals with the 
estimation of the adsorption properties of Amberlite FPA-98 that constitutes a possible source of adsorbent that 
could be used for the removal of dyes from textile wastewater and, more generally, in industrial wastewater . 
The focus of the present study was to assess the potentiality of Amberlite FPA-98 as an adsorbent for the dye 
A07 from aqueous solution as an ideal alternative to the current expensive methods of removing dyes from 
wastewater using the treatment of a synthetic textile effluent containing an acide dye, Acide Orange 7 (A07), as 
the adsorbate. The Amberlite FPA-98 will be, first, characterized in terms of chemical composition, structure 
and texture [10]. 

Adsorption studies were carried out under various parameters such as flow rate, initial dye 
concentration, and bed height. The continuous adsorption in fixed-bed column is often desired from industrial 
point of view. It is simple to operate and can be scaled-up from a laboratory process [11]. A continuous packed 
bed adsorber does not run under equilibrium conditions and the effect of flow condition (hydrodynamics) at any 
cross-section in the column affects the flow behaviour downstream. Breakthrough determines bed height and the 
operating life span of the bed and regeneration times [12]. Adsorption in fixed-bed columns using activated 
carbon has been widely used in industrial processes for the removal of contaminants from aqueous textile 
industry effluents, since it does not require the addition of chemical compounds in the separation process [13]. 
Adsorption in a fixed bed column can be used continuously under high effluent flow rates and it has been used 
in many pollution control processes such as removal of ions by an ion exchange bed or removal of toxic organic 
compounds by carbon adsorption [14]. 

In this study, the amberlite has been tested for removal of aqueous solutions. In this paper, 
amberlite is used to remove azo dye (A07) from aqueous solution through column studies. The objective of this 
study was to investigate the adsorption potential of A07 onto amberlite fixed-bed. The important design 
parameters such as inlet concentration of dye solution, fluid flow rate and column bed height [15] were 
investigated using a laboratory scale fixed-bed column. 

In application of adsorption process on an industrial scale, it is crucial to improve process efficiency, 
reduce operational cost and time to minimum and take into account the most important factors, what can be 
achieved by applying the optimization techniques such as response surface methodology (RSM). Determining 
the effect of a single factor on the efficiency of the process is relatively simple . It is definitely more of a 
challenge to assess the effect of several parameters at once. Response surface methodology based on 
experimental data makes easier to plan the entire modeling process by reducing the number of experiments to 
the necessary minimum, and allows a mathematical equation to fit the experimental results, which is required 
for the process optimization [16].In general, RSM is a mathematical technique applied in the progression of an 
appropriate functional relationship between the response and the related input variables. The structure of this 
relationship is unknown, but can be approximated by low-order polynomials (the most common are first and 
second-degree polynomial models) . This methodology helps to determine the most important parameters and its 
main, quadratic effect or interactions which influence the response. RSM has been extensively used as an 
optimization, prediction and interpretation technique for factorial designs [17].RSM is an useful tool for the 
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modeling and analysis of systems in which response of interest depends on several factors and the relationship 
between independent and dependent variables in a system is unknown. RSM modeling has been successfully 
applied for biosorption in the past few years . RSM was selected as an effective statistical and mathematical 
approach in order to recognize the efficiency of an experimental system [18]. Various parameters were 
simultaneously appraised using RSM with a minimum number of experiments. Therefore, a study conducted by 
RSM can reduce the cost, decrease process variability and need less time in comparison to the conventional one- 
factor-at-a-time statistical strategy [19]. 

The present study investigates the application of response surface methodology approaches to 
predict adsorption capacity of amberlite for the removal of A07 from aqueous solution. The effect of various 
operational parameters, including the initial dye concentration, flow rate, and bed height on the adsorption of 
A07 were examined and modeled by response surface methodology RSM. The Box-Behnken design BBD and 
the full factorial design FFD with structure optimized by RSM models were compared in terms of predictability 
and accuracy of fit, taking into account their implementations and limitations. Optimization of BBD by response 
surface method is completely novel approach of RSM approximation application in chemical processes. The 
main objectives of this work are to investigate the individual and the interactive effects of three operating 
parameters, mainly: initial dye concentration, flow rate, and bed height on the adsorption of A07 in a fixed bed 
reactor by using a BBD and full factorial design (FFD) [20], 

Conclusions as a result, the adsorptions of dyes onto AC were commonly investigated using 
traditional methods, but the A07 adsorptions based onto amberlite FPA 98 have not been studied. Moreover, 
Box-Behnken design (BBD) is rarely used for the A07 removal from aqueous solution. Then in the text, the 
adsorptions of dye onto amberlite FPA 98 have been investigated in continuous systems. The main effects and 
interaction effects between process variables on the dye adsorptions were analyzed based on the BBD. Their 
maximum adsorption capacities have been optimized using RSM method. [21], 

II. Material and Methods 

A. Porous Polymers (Resin Adsorbent) .Amberlite FPA 98 

There are two main criteria for the choice of adsorbent: The maximum adsorption capacity (q IT1 ): 
solute concentration transmitted in the solid phase, must be maximized. The adsorption kinetics: the adsorption 
must take place rapidly. Exchange resins, strongly basic anionic ion (85% adsorbent and 15% anion exchange), 
the principle is to exchange certain ions, or all with active groups on the resins. Resins, which are in the form of 
beads or powder form are thus able to exchange mobile ions with ions of the same sign, contained in a solution 
with which they are contacted. Strongly basic anionic resins bind the strong and weak acid anions in a pFl range 
between 1 and 12. This type of resin is regenerated with sodium hydroxide solution and in the form of grains 
aspheric 0.5 to 1.5 mm in diameter, with specific surface areas of up to 750 m 2 /g. 
Amberlite FPA98 is unique, food-grade acrylic macroporous strong base anion resin. This use of ion exchange 
technology based discoloration was more effective and more economical. The application of these adsorbents 
AMBERLITE FPA98 Cl mainly develop in water treatment in the domain refining of sugar as adsorbents dyes 
capturing and purification of pharmaceuticals and The guard beds for precious chromatography media [22]. 

B. Adsorbate Preparation 

The dye chosen in this study is the Acide Orange 7, also called Acide Orange II (Sigma-Aldrich), 
belonging to the family of the anionic dyes. It is representative of a textile type of pollution. Its main features 
are represented in Table I; its structural formula is shown in Figure 1. Stock solutions were prepared by 
dissolving requisite quantity of dye without further purification in distilled water, and the concentrations used 
were obtained by dilution of the stock solution. The pH was adjusted to a given value by addition of HC1 (IN) 
or NaOH (IN) [23], 



Fig. 1 Molecular structure of Acide Orange 7 
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TABLE I: MAIN CHARACTERISTICS OF THE BASIC DYE ACIDE ORANGE II 


Dye 

Acide Orange 7 

Molecular formula 

C16H11 N2Na04S 

Molecular weight [g/mol] 

350.33 

Molecular volume (A/molecule) 

231.95 

Molecular surface area (A/molecule) 

279.02 

Width (A) 

7.3 

Length (A) 

13.6 

Depth (A) 

2.3 

X (nm) 

485 

pKa 

pKl 11.4; pK2 1.0 


III. SORPTION EXPERIMENTS IN FIXED-BED TECHNIQUE 

A. Experimental Procedures 

The fixed bed experiments were carried out in a glass column of 2 cm internal diameter, 30 cm of the length 
height and three sampling points at 5 cm intervals. A known quantity of Amberlite FPA-98 was packed in the 
column to yield the desired bed height of the adsorbent 50, 100 and 150 mm (equivalent to 3.5, 7 and 10,5 g of 
activated carbon) with a layer of glass wool at the bottom. Distilled water was passed through the column in 
order to remove impurities from the adsorbent. Three flow rates (2, 4 and 6 mL/min) were pumped to the top of 
the packed column by using peristaltic pump with different initial dye concentrations (5, 30, 50, 80 mg/L). The 
samples of A07 solutions at the outlet of the column were taken at regular time intervals and the concentration 
of dye was measured using an UV-visible spectrophotometer (Neptune Chemical Pump) at wavelength of 485 
nm. Fixed bed studies were terminated when the column reached exhaustion [24]. 

The schematic diagram of fixed bed column used in adsorption study is shown in Fig. 2. The experimental detail 
is given in Table II. Briefly, the experiment was carried out by passing through column (packed with 21 g of 
Amberlite FPA-98) with controlled flow rate . 


I l 

-*-l 

i I 



|_| Feed Storage 

Fig. 2 Schematic diagram of fixed bed column used in adsorption study of A07 onto Amberlite FPA 98. 


TABLE II 

EXPERIMENTAL DETAILS FOR COLUMN ADSORPTION OF A07 ONTO AMB FPA 98 


EFFECT OF 




SYSTEM 

Flow rate(mL/min) 

[A07](ppm) 

Bed height(cm) 

Flow rate 

2, 4 and 6 

30 

15 

Initial Concentration 

15 

10, 30 and 50 

15 

Bed height 

15 

30 

5, 10 and 15 
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B. Removal efficiency 


The dye removal percentage using Amberlite FPA98 adsorbent was calculated from: 


Re(%) = lOO x (C 0 - C = )/C Q 


(1) 


Where Co is the initial concentration of dye in solution (mg/L), and C e is the final dye concentration in aqueous 
solution after phase separation (mg/L). 

C. Response surface methodology 

Response surface methodology is an experimental technique used for predicting and modeling 
complicated relation between independent factors and one or more responses. [25].In this study, response 
surface methodology was applied to optimize the adsorption of Acide orange 7 by AmberliteFPA98 . 
Experiments were performed using Box- Behnken design (BBD) and full factorial design (FFD). The second- 
order polynomial equation extended with additional cubic effects was employed as an objective function. The 
second order model is usually sufficient for the modeling and optimization on the basis of designs, however 
third order and higher effects are sometimes important, especially in order to achieve better fit and insignificant 
lack-of-fit. For instance, Box-Behnken design was created to estimate the second-order model, however there 
may be situations in which non- random portion of this model provides an inadequate representation of the true 
mean response, an indication of lack-of-fit of the second- order model. Thus, in this study some third order 
model terms were added to the second order polynomial equation. Accuracy of model fitting was evaluated by 
means of ANOVA. All calculations were performed in Statgraphics Software [26]. 

I j .I - -, 

D. Box-Behnken design 

1 ■ _ ■“‘ tl | ■ * _ —1 

In this study, the BBD design methodology was employed to optimize the operational variables and 
was used to predict impacts of respective parameters on the adsorption process. Among many factors affecting 
the adsorption process, three process variables , i.e. initial A07 dye concentration (Xi), bed height (X 2 ) and flow 
rate (X 3 ) were selected and were considered as independent variables and the removal of dye (Y) as a response 
was defined and modeled. BBD contains set of 15 experiment runs whose values of each factor with three 
levels (low, medium ,high), being is coded to standard values (-1, 0,+ l) in the appropriate range and levels of 
parameters were listed in (Table III). The parameters ( X,) were coded as x, via Eq. (2): 


( 2 ) 


x = (x, -x Q ) 

1 _srv" 


where X 0 and dX are the values of X, at the center point and step change, respectively . The second-order 
polynomial response equation was used to probe the interaction between the dependent and independent 
variables.The removal (%)of dye is selected as the response for the combination of independent variables. 
Subsequently experimental data was fitted to the second order polynomial model extended with additional cubic 
interaction effects (Eq. (3)) using the least square procedure as follows: 

Y = b 0 + b,x, + b 2 x 2 + b 3 x 3 + b,,x ; + b l2 rfx 2 + b x 3 v ( x 3 + b 22 x 2 + b 23 x 2 x 3 + b 33 x 3 + e (3) 

n i |Y Vi i j f' L ~v 

where Y is the predicted response associated with each factor level combination ;The coefficients in the 
equation represent: the intercept (b 0 ) is constant, the main (bi,b 2 , b 3 ) are linear effect, quadratic effect (bn, b 22 , 
b 33 ) and interactions (bi 2 , b 23 , b ]3 ) effects, respectively; x and Xj are the coded values of independent variables; 
and e is the residual error. Validation of the model fit and significance analysis of variables were performed 
using analysis of variance (ANOVA). The results were analyzed by analysis of variance (ANOVA) and a 
calculation correlation coefficient (R 2 ) between predicted and experimental points [27], 

Nonlinear X 2 analysis is a useful method, which can compare the experimental and model predicted data. And it 
is estimated using the following equation: 



(4) 


whereY exp . (%) and Y pred . (%) are the adsorption removal of A07 onto Amberlite FPA 98 by experiment 
determined and model predicted. 
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E. Experimental design 

To determine the optimal conditions for the main parameters, a Box-Behnkendesign (BBD) was applied. For 
the adsorption process, significant variables, such as the initial dye concentration, flow rate, and bed height were 
regarded as the independent variables and designated as Xj-X 3 , respectively. The dye concentration (X|j range 
of 10-50 mg/ L, flow rate (X 2 ) range of 2-6 ml/min, bed height (X 3 ) of 5-20 cm were chosen as given in Table 
III. 


TABLE III : INDEPENDENT PROCESS VARIABLES AND THEIR EXPERIMENTAL LEVELS USED FOR BOX-BEHNEKEN 

DESIGN (BBD). 


Variables, unit 

Factors 


Levels 



X 


-i 

0 

+ 1 

Initial dye concentrations of A07 (mg/L) 

X, 


10 

30 

50 

bed height (cm) 

X 2 

r — _ 

5 

10 

15 

flow rate (ml/min) 

X 3 

~ L '■ 

2 

4 

6 


F. Selection of the significant parameters 


The Box-Behnken design consists of 15 experimental points. The experimental conditions, and the 
adsorption capacity obtained for each point set by the Box- Behnken design are shown in Table IV (1-11), 
together with the three central point repetitions (12-15). The relationship between responses and processed 
variables was examined for the response approximation function (Y) using Eq. (1), following by the statistical 
analysis of the model obtained. The most significant process variables were identified by Box-Behnken design 
(BBD) experimental design. The advantage of this design is its ability to investigate of a large number of 
factors in a relatively low number of experimental runs, hi this study 15 run BB design with 3 factors, including 
A07 dye concentration (Xj), flow rate (X 2 ) and bed height (X 3 ) was considered. Each independent variable was 
tested at two levels, high and low, which were -1 and +1, respectively. All experiments were conducted in 
duplicate and the average values of adsorption capacity were taken as a response (Y). The matrix design is 
shown in Table IV. On the basis of BBD three the most significant parameters were chosen for further 
investigation (modeling and optimization by RSM ). 

TABLE IV : BOX-BEHNKEN DESIGN MATRIX WITH CODED AND UNCODED VALUES OF THE INDEPENDENT VARIABLES 
INFLUENCING ADSORPTION OF A07 BY AMBERLITE WITH EXPERIMENTAL AND PREDICTED VALUES OF THE 

RESPONSE. 


Experimental run No. 

Coded values (uncoded values) 

Removal dye Y [(%)] 

X] [ppm] 

X 2 [cm] 

X 3 [ml/min] Experimental values 

Predicted responses 

i 

0(30) 

-1 (5) 

-1(2) 

50 



53,2375 

2 

0(30) 

1 (5) 

1(6) 

66 



64,7375 

3 

1 (50) 

0(10) 

1(6) 

69 

J 1 


70,0125 

4 

-1 GO) 

-1 (5) 

0(4) 

"48,9 



46,675 

5 

-1 GO) 

0(10) 

-1(2) 

52 



50,9875 

6 

0 (30). . .. 

f (15) 

-1(2) 




59,2625 

7 

1 (50) \ 

V (15) 

-1 (2) 

1 



34,5125 

8 

0(30) 

1 (15) 

1(6) 

79 



75,7625 

9 

'-l GO) 

0(10) 

1(6) 

40 



43,4875 

10 

1(50) 

1 (15) 

0(2) 

58 



60,225 

11 

1 (50) 

-1 (5) 

0(2) 

49 



49,25 

12 

-1 GO) 

1 (15) 

0(2) 

53 



52,75 

13 

0(30) 

0(10) 

0(2) 

58 



57 

14 

0(30) 

0(10) 

0(2) 

57 



57 

15 

0(30) 

0(10) 

0 (2) 

56 



57 


G. Full factorial design 


Factorial designs allow the simultaneous study of the effects that several factors may have on the 
optimization of a particular process. It determines which factors have the important effects on the response as 
well as how the effect of one factor varies with the level of the other factors. The effects are the differential 
quantities expressing how a response changes as the levels of one or more factors are changed. Also, factorial 
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designs allow measuring the interaction between each different group of factors. The interactions are the driving 
force in many optimizations of the processes. Without the use of factorial experiments, some important 
interactions may remain undetected, and the overall optimization may not be attained. One of the simplest types 
of factorial designs used in experimental work is one having two levels (2 k ). In a 2 k factorial design experiment, 
each factor may be assigned two levels: low (-1) and high (+1). If k factors are considered, then 2 k 
measurements are required to perform a factorial design analysis [28]. In this investigation, three operating 
factors were chosen as independent variables, namely: initial A07 dye concentration (Xi), flow rate (X 2 ) and 
bed height (X 3 ) was used to predict the removal of dye . The natural values of each factor and their respective 
levels are presented in Table V. 


TABLE V : THE EXPERIMENTAL RANGES AND LEVELS OF INDEPENDENT VARIABLES. 


Levels 

[A07]o (mg/L) 

Qv (mL/min) 

H (cm) 



r" 


-i 

10 

1 2 

.J 1 1. I 

5 

+i 

50 

L 

20 


The selection of levels of different factors is carried out on the basis of the preliminary trials and previous 
publishing result: Initial dye concentration [A07]o ranging from 10 to 50 mg/L, flow rate from 2 to 6 mL/min 
and bed height from 5 to 20 cm . The design performed according to Table VI was composed of 2 3 factorial 
designs. 




3 1 



TABLE VI: 2 3 FULL FACTORIAL DESIGN. FULL FACTORIAL DESIGN WITH CODED AND UNCODED VALUES OF THE 
INDEPENDENT VARIABLES AND EXPERIMENTAL AND PREDICTED VALUES OF THE RESPONSE. 




Coded values (uncoded values) 

Removal dye Y [(%)] 

run No. 

X, [ppm] 

X 2 [ml/min ] 

X 3 (cm) Experimental values 

Predicted 

responses 

i 

-1 (10) 

-1 (2) 

-1(5) 

55 


54,625 

2 

1(50) 

-1 (2) 

1 (20) 

57 


56,625 

3 

-1 GO) 

-1 (2) 

1 (20) 

52 


52,375 

4 

1(50) 

1(6) 

-1 (5) 

56 


55,625 

5 

1(50) 

1(6) 

1(20) 

50 


50,375 

6 

1(50) 

-1 (2) 

-1 (5) 

62 


62,375 

7 

-1 (10) 

1 (6) 

1(20) 

53 


52,625 








8 

-1(10) 

—fr(6) 

-1 (5) 

54 


54,375 


The experimental data are analyzed by full factorial design to fit the following first order polynomial equation 
(Eq.(5)). 


Y — b 0 + bjXj -T b 2 x 2 -+- b,x 3 -T h ] -t- h ] . A, x . -t- b 23 x 2 x 3 -t- b l23 x 1 x 2 x 3 


+ e 


(5) 


Where Y the estimated removal of dye ; b 0 is the value of fitted response at the centre point of design; bj and bjj 
are the linear and interaction terms, respectively [29].When the response data are obtained from the test work, a 
regression analysis is carried out to determine the coefficients of the response model (bi; b 2 ; . .. b n ), as well as 
their standard errors and their significance. In addition to the constant (b 0 ) and error (e) terms, the response 
model incorporates .Linear terms in each of the variables (xi; x 2 ; . . . ; x n ). First-order interaction terms for each 
paired combination (xix 2 ; XiX 3 ; . . . ; x n _iX n ) 


H. Analysis of variance (ANOVA) 
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ANOVA expounds every variation in the statistically obtained model and importance of each model 
parameters. The significance of the model was evaluated by F-test for a confidence level of 95% as well as lack- 
of-fit test. In general, the greater the F-value and the smaller the p-value, the more significant is a model. 
Moreover, effects and their importance in the model were investigated adapting t-test and p-value. Usually, the 
larger the t-value and lower probability p-value (p < 95%), the model parameter is considered as significant 
[30]. The sum of squares, degree of freedom and mean squares were also determined for the model and error. 


IV. Results and discussion 


A. Box-Behnken 

Response surface methodology (RSM) is more advantageous than the traditional single parameter 
optimization because it can save time, space and raw material. In experimental design, a Box- Behnken design 
(BBD) is a type of RSM, and it is used for optimizing the important process variables. The most important 
parameters, which affect the efficiency of adsorption of A07 onto AMB FPA 98 , are A07 dye concentration, 
flow rate and bed height of the solution in a continuous fixed bed. In order to study the combined effect of these 
factors, experiments are performed for different combinations of the physical parameters using statistically 
designed experiments. The initial dye concentration range studied is between 10 and 50 ppm. The flow rate is 
between 2 and 6 ml/min. The bed height is varied between 5 and 20 cm . The main effects of each of the 
parameter on A07 removal efficiency is given in Figs.3. Fig. 3 shows that the removal efficiency increases 
with increasing A07 dye concentration, flow rate or bed height. Consequently, we note high A07 removal 
efficiencies at high flow rate and bed height. 


Main Effects Plot for removal dye 



- 1.0 1.0 - 1.0 1.0 - 1.0 1.0 
intial concentration bed height flow rate 

‘ I I '■ L 

Fig. 3. Main effects plot of parameters for A07 removal efficiency. 

Table VII shows the experimental results of removal of A07 in the solution for the 15 experiments. Using the 
experimental results, the regression model equations (second-order polynomial) relating the response is 
developed and is given in Eqs. (6). Apart from the linear effect of the parameter for the response, the RSM also 
gives an insight into the quadratic and interaction effect of the parameters. These analyses are done by means of 
Fisher’s ‘F’ test and Student ‘t’ test. The Fisher’s ‘F’ test is used to determine the significance of each of the 
interaction among the variables, which in 

turns may indicate the patterns of the interactions among the variables. 

'. I i L i f * I_ r j' l " — ? 

TABLE VII THE EXPERIMENTAL DATA FOR A07 REMOVAL EFFICIENCY IN SOLUTION ACCORDING TO BBDESIGN. 


\ 

a v ,i i r \ 



ci-T-b ... / 

Removal dye Y [(%)] 


Experimental run No. 

Experimental values 

Predicted responses 


i 

50 


53,2375 

2 

66 


64,7375 

3 

69 


70,0125 

4 

48,9 


46,675 

5 

52 


50,9875 

6 

58 


59,2625 

7 

38 


34,5125 

8 

79 


75,7625 

9 

40 


43,4875 

10 

58 


60,225 

11 

49 


49,25 
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12 

53 

52,75 

13 

58 

57 

14 

57 

57 

15 

56 

57 


In general, the larger the magnitude of F, the smaller the value of P, the more significant is the corresponding 
coefficient term. The regression coefficient, F and P values for all the linear, quadratic and interaction effects of 
the parameter are given in Tables VIII for the removal of A07. It is observed that the coefficients for the linear 
effect of the factors flow rate and bed height (P =0.0191and p=0.0025) for the responses are significant except 
A07 dye concentration (P = 0.1008) for removal of dye is slightly less significant. 

However, for the removal efficiency of A07, the interaction effect of the variables dye concentration and flow 
rate is found highly significant p=0.0017 exempt the interaction between dye concentration and bed height (P = 
0.5194). Consequently, the best fitting response function , for the A07 removal efficiency model are then 
conveniently written as follows: 

Y = 57,0 + 2,5125.X! + 4,2625.x 2 + 7,0.x 3 - 9,1375.x, 2 + l,225.x,.x 2 + 10,75.x,.x 3 + 4,3625.x 2 2 + l,25.x 2 .x 3 + 
1,8875.x 3 2 (6) 

l' "i 1 i 

where Y (%) is the removal dye of A07 , and x,, x 2 and x 3 are the A07 dye concentration , flow rate and bed 
height respectively. 

TABLE VIII ANALYSIS OF VARIANCE AND CORRESPONDING F AND P VALUES FOR A07 REMOVAL EFFICIENCY. 


Source 

Sum of Squares 

Df 

Mean Square F-Ratio 

P-Value 

A:intial concentration 

-50,5012 

1 

50,5012 

4,04 

0,1008 

B:bed height-* 

145,351 

1 

145,351 

11,62 

0,0191 

C:flow rate 

392,0 

1 

392,0 

31,33 

0,0025 

AA 

308,285 

1 

308,285 

24,64 

0,0042 

AB 

6,0025 

1 

6,0025 

0,48 

0,5194 

AC 

462,25 

1 

462,25 

36,95 

0,0017 

BB 

70,2698 

1 

70,2698 

5.62 

0,0640 

BC 

6,25 

1 

6,25 

0,50 

0,5113 

CC 

13,1544 

‘vi 

13,1544 

1,05 

0,3522 

Total error 

62,5525 

5 

12,5105 



Total (corr.) 

1547,04 

14 





The ANOVA table partitions the variability in removal dye into separate pieces for each of the effects. It 
then tests the statistical significance of each effect by comparing the mean square against an estimate of the 
experimental error. In this case, 4 effects have P-values less than 0,05, indicating that they are significantly 
different from zero at the 95,0% confidence level. The R-Squared statistic indicates that the model as fitted 
explains 95,9566% of the variability in removal dye. The adjusted R-squared statistic, which is more suitable 
for comparing models with different numbers of independent variables, is 88,6785%. Further, the ANOVA for 
A07 removal efficiency in solution indicates that the second-order polynomial model Eqs. (6) is highly 
significant and adequate to represent the actual relationship between the response and variables, with very a 
high value of coefficient of determination (R 2 = 0.959566 for A07 removal efficiency in solution. This implies 
that 95.95% of sample variation for A07 removal efficiency in solution is explained by the model. 

The statistical significance of the ratio of mean square variation due to regression and mean square 
residual error is tested using ANOVA. ANOVA is a statistical technique that subdivides the total variation in a 
set of data into component parts associated with specific sources of variation for the model. According to the 
ANOVA (Table VIII), the F statistics values for all regression are higher. The large value of F indicates that most of 
the variation in the response can be explained by the regression equation. The associated P value 
is used to estimate whether F stat i st i cs is large enough to indicate statistical significance. 

The ANOVA table also shows a term for residual error, which measures the amount of variation in the response 
data left unexplained by the model. The form of the model chosen to explain the relationship between the factors 
and the response is correct. 
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The 3D response surface and 2D contour plot are generally the graphical representation of the regression 
equation. We will use it to search the optimal values of the process parameters. Then, the response surface plots 
and contour plots to estimate the removal efficiency (Figs. 4 and 5) is given. 


t 

1 



Fig. 4. Response surface plot of A07 removal efficiency . 






Contours of Estimated Response Surface 
flowrate=0,0 



Fig. 5. Contour plot of estimated response surface of A07 removal efficiency. 


Thus, the surface and contour plots for A07 removal efficiency in Fig. 7 shows the interaction 

effect of bed height and initial concentration .The response surface of mutual interactions between the variables 

is found to be elliptical and the maximum A07 removal efficiency is obtained in the following cases: 


The bed height and initial concentration increase simultaneously. 

The initial concentration increases and the bed height is between 10 and 15 cm and 
remains unchanged. 

The bed height increases and initial concentration is between 10 and 20 PPM and remains 
stable. 

■- i. _ / } 

The geometrical representation of the response removal dye , when The bed height and initial concentration 
increases the removal dye increases. We also note that, the influence of flow rate is not significant. Then, to 
have a good removal dye it is beneficial to work with high bed height of column . The highest value of the bed 
height which gives maximum of A07 removal is 15 cm. 

\ h? I. 


B. Full factorial design 


i 


The model equation for adsorption of A07 by amberlite in fixed bed was obtained after performing eight 
experiments and discarding the insignificant effect (bl2) is as follows using some statistical tests . Table IX 
shows the experimental results of removal efficiency of A07 in the solution for the 8 experiments. 


TABLE IX THE EXPERIMENTAL DATA FOR A07 REMOVAL EFFICIENCY IN SOLUTION ACCORDING TO FF DESIGN. 


run No. 

Removal dye Y [(%)] 

Experimental values 

Predicted responses 

i 

55 

54,625 

2 

57 

56,625 

3 

52 

52,375 

4 

56 

55,625 
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5 

50 

50,375 

6 

62 

62,375 

7 

53 

52,625 

8 

54 

54,375 


Y=54,875 + l,375.xi - l,625.x 2 - l,875.x 3 - 1,625.xi.x 2 - 0,875.xi.x 3 + 0,125.x 2 .x 3 


(7) 


The model’s coefficients were estimated by standard lest square regression techniques using an EXCEL 
software. A good adjustment of the (Eq. (7)) to the experimental data was verified through the high correlation 
coefficient value obtained R2 = 0.987 (Fig. 2). The (Fig. 2) shows the absence of a trend, indicating that the 
mathematical model is adequate and that there is no inconsistency between the experimental and calculated 
values of the response. 






Vi 


r /' 


w 



- -c_ 


Fig.6. Analysis of quality of model by comparison of experimental and predicted responses. 


fig.6 shows that the difference between the measured and the predicted values do not exceed 1%. Therefore, 
all those results indicate that the model can adequately represent the data. Initial A07 concentration [AO7] 0 (xi) 
has the strongest effect on the response since coefficient of x 3 (bl 6= +1,375) is large than the coefficients of the 
other investigated factors, the positive sign indicate that there is a direct relation between Initial 
A07concenU'ation [AO7] 0 and response . According the regression equation, flow rate (x 2 ) and bed height (x 3 ) 
have a negative effect on the response (b 2 = -1,625 and b 3 = - 1,875) respectively which have been explained by 
reduction in residence time. The significance interactions found by the design of experiments are between flow 
rate-bed height (x 2 x 3 ) with effect (b 23 = +0,125). 

\ j,i ! 1 C f lL ~ £ 

The main effects of each of the parameter on A07 removal efficiency is given in Figs.7. Fig. 7 shows that 
the removal efficiency increases with increasing A07 dye concentration and decreases with increasing flow rate 
or bed height. Consequently, we note high A07 removal efficiencies at high A07 dye concentration and low the 
both of flow rate and bed height 
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Main Effects Plot for removal dye 



Fig. 7. Main effects plot of parameters for A07 removal efficiency. 

Results of analysis of variance were summarized in Table X. They indicated that the three regressions, F and P values for the linear and 
interaction effects of the parameter were significant at the probability level of 95 % to 99 %. 


■S, L — --* i 

TABLE X ANALYSIS OF VARIANCE FOR REMOVAL DYE 


Source 

Sum of 
Squares 

Df 

Mean 

Square 

F-Ratio 

P- 

Value 

A: initial dye 

15,125 

1 

15,125 

13,44 

0,1695 

concentration 






B:bed height 

21,125 

1 

21,125 

18,78 

0,1444 

C:flow rate 

28,125 

1 

28,125 

25,00 

0,1257 

AB 

21,125 

1 

21.125 

18,78 

0,1444 

AC 

6,125 

1 

6,125 

5,44 

0,2578 

BC 

0,125 

1 

0,125 

0,11 

0,7952 

Total error 

1,125 

1 

1,125 


\K 

Total (corr.) 

92,875 

7 


/ 



The ANOVA table partitions the variability in removal dye into separate pieces for each of the effects. It then 
tests the statistical significance of each effect by comparing the mean square against an estimate of the 
experimental error. In this case, 0 effects have P-values less than 0,05, indicating that they are significantly 
different from zero at the 95,0% confidence level. The R-Squared statistic indicates that the model as fitted 
explains 98,7887% of the variability in removal dye. The adjusted R-squared statistic, which is more suitable 
for comparing models with different numbers of independent variables, is 91,5209%. 

,J I * \ I 'll 

The 3D response surface and 2D contour plot are generally the graphical representation of the regression 
equation. We will use it to search the optimal values of the process parameters. Then, the response surface plots 
and contour plots to estimate the removal dye of A07 (Figs. 8 and 9) is given. 
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Estimated Response Surface 
flow rate=0,0 



Fig. 8. Response surface plot of A07 removal efficiency . 


\ \ 
—iJl 


Contours of Estimated Response Surface 
flow rate=0,0 



initial dye concentration 


removal dye 
— 50,0 


- 53,0 

- 54,5 
56,0 


Fig. 9. Contour plot of estimated response surface of A07 removal efficiency. 


Thus, the surface and contour plots for A07 removal efficiency in Fig. 9 shows the interaction 
effect of bed height and initial dye concentration .The response surface of mutual interactions between the 
variables is found to be elliptical and the maximum A07 removal efficiency is obtained in the cases: The initial 
concentration increase and the bed height decrease simultaneously. The initial concentration increases and the 
BED F1EIGHT is between 5 and 10 cm and remains unchanged. 

The bed height decreases and initial concentration is between 40 and 50 ppm and remains stable. 

The geometrical representation of the response removal dye, when The bed height decreases and initial 
concentration increases the removal dye increases. We also note that, the influence of flow rate is not 
significant. Then, to have a good removal dye it is beneficial to work with low bed height of column. The 
lowest value of the bed height which gives maximum of A07 removal is 5cm. 


V. Conclusion 




The present study clearly demonstrated the applicability of adsorption process using the fixed bed for A07 
removal. This study clearly showed that RSM is one of the suitable methods to optimize the best operating 
conditions to maximize the A07 removal. BB and FF design is successfully employed for experimental design 
and analysis of results. The Amberlite FPA 98 , which was used without further purification for the removals 
of A07 from aqueous solution because it leans close to practical purposes. The process variables of removal of 
dye by amberlite FPA 98 have been optimized based on RSM method and the individual and interaction effects 
of the process variables were investigated. The results indicated that all the three process variables have a direct 
relationship for the removal A07 dye onto amberlite. Satisfactory empirical model equations are developed 
for the removal of A07 in solution using RSM to optimize the parameters. Graphical response surface and 
contour plot is used to locate the optimum point. 
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